1 



cn 

> 
o 



Geometric Measurement of Topological Susceptibility on Large LatticesB 

J. Grandy and R. Gupta*^ 

'"T-S Group, MS B285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 U. S. A. 

The topological susceptibility of the quenched QCD vacuum is measured on large lattices for three (3 values 
from 6.0 to 6.4. Charges possibly induced by 0(a) dislocations are identified and shown to have little effect on 
the measured susceptibility. As /3 increases, fewer such questionable charges are found. Scaling is checked by 
examining the ratios of the susceptibility to three physical observables computed in Refs. [|l|, ^, ^, ^ and Aiatt- 
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The topological susceptibility of the QCD vac- 
uum is defined by 
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d'^xF{x)F{x) . (1) 



The susceptibility is related to the -q' and rj masses 
in the chiral large- A^c limit through the Witten- 
Veneziano formula [ ^ 
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The computation of xt on a lattice is not 
straightforward. The simplest method, replacing 
FF by a product of plaquettes, suffers mixing 
with F'^ and a constant [ |^, and a large pertur- 
batively determined renormalization factor, but 
cooling and smearing algorithms have been pro- 
posed to refine this method [ |^, ^. The geomet- 
ric method, first proposed by Liischer[ ^, is not 
susceptible to multiplicative renormalization but 
contains additive divergences from compact dis- 
locations [ |l^ so that 
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The dislocation described by Gockeler et aL[ 
produces a power po < in (^) and thus a power 
law divergence in the measured susceptibility as 
a ^ 0. The factor C{p) for p < is crucial to de- 
termining the onset of divergent behavior in the 
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measured susceptibility for finite a, and we inves- 
tigate whether divergent terms affect the result 
at realistic values of a used in our calculation. 

We are examining several aspects of topology 
on the lattice. First, we compute xt using the geo- 
metric method and examine the systematic effect 
of dislocations in the range /? = 6.0 to (3 = 6.4. 
We plan to study the various methods of mea- 
suring Xf, including cooling, and use the geomet- 
ric topological charge density to examine the spin 
asymmetry in deep inelastic scattering off a pro- 
ton. At present we have completed the first stage 
and report the results here. 

Our lattice ensembles are listed in Table 0. We 
evaluate the integrals derived by Gockeler et al.[ 
on the surfaces of each hypercube on the lat- 
tice, obtaining a local integer for each hypercube, 
and add the local integers in each configuration 
to obtain the topological charge Q. The integra- 
tion requires a two-stage process. First, all hy- 
percubes are integrated using two different coarse 
(~ 5^) cubic meshes. Hypercubes which show no 
definite evidence of convergence toward zero are 
integrated more accurately, with ~ 10^ meshes on 
subcubes. Most hypercubes are brought within 
10^^ to 10~^ of an integer with the finer mesh, 
but a few difficult hypercubes require further re- 
finement. The number of difficult hypercubes per 
unit topological charge decreases as /3 increases. 

We attempt to separate spurious charges, 
caused by low-action dislocations, from physical 
charges. These dislocations in SU{2) typically 
contain a central plaquette near —1 and have a 
radius of about 1 lattice spacing [ Since the 
topological structure of SU{3) is determined by 
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Figure 1. Normalized nearest neighbor correla- 
tions from charged hypercubes to nearest plaque- 
tte with RDC < 0.15 (Eq. |), versus r^,„, in lat- 
tice units. Highlighted correlations at rmm < 2a 
indicate possible compact dislocations. 



embedded SU (2) windings, we search for correla- 
tions between hypercubes of nonzero charge and 
plaquettes near the SU{3) cut locus, the analog 
to the —1 group element in SU{2). Plaquettes 
whose radial distance to cut locus (RDC)[ [l^, 

RDC=(^2g.?j (_--) (4) 

is less than 0.15 are selected for presentation here. 
In this definition (^) the di are the three eigenan- 
gles of the SU{3) matrix in ascending order, so 
that their sum is zero. For each charged hyper- 
cube in the ensemble the distance rmin from the 
center of the hypercube to the center of the near- 
est selected plaquette is found. The number of oc- 
currences at each lattice distance rmin, N{r,nin) 
is divided by the weight factor g{rmin), the num- 
ber of available plaquettes at distance rmin, to 
obtain a normalized nearest neighbor correlation 



Figure 2. Xt^VV^ (circles), xyViWOAiatt) 
(crosses), xj^/vip (diamonds), and Xt^''/(10/7r) 
(squares) versus /3. Scaling is indicated when the 
ratios are independent of (3. The actual values of 
these ratios are unimportant, since we consider 
only their dependence on /3. The IG'^ x 40 lattice 
results are shown for /3 6.0. 



n{rmin)- If the charges and nearly critical pla- 
quettes are decoupled, n{rmin) will fall exponen- 
tially as a function of volume spanned by rmin, 

As we investigate possible dislocations, we find 
some encouraging results. Our plot (Fig. |l|) 
shows n(rmin) decaying exponentially as a func- 
tion of r^in , except for a dramatic rise for rmin < 
2a, indicating that questionable charges within 
2a of a nearly critical plaquette may indeed be 
dominated by compact dislocations which are lat- 
tice artifacts. Where Nq is the number of ques- 
tionable charges in a configuration, we find that 
{Ng)/{Q'^) = 0.21 at ^ = 6.0, and 0.05 at ^ = 6.4. 
This decrease suggests that for the /? values con- 
sidered here the additive term in Eq. ^ is dom- 
inated by the p > part of the integral, and 
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Table 1 



Lattice Ensembles and Results 





Lattice 


Sample 


xr(v^) 


Xt (wp) 




Size 


Size 


(MeV) 


(MeV) 


6.0 


16^ X 40 


34 


228(15) 


264(19) 


6.0 


24^ X 40 


23 


220(29) 


255(34) 


6.2 


32^ X 48 


22 


251(16) 


251(28) 


6.4 


323 X 48 


21 


268(18) 


264(20) 



the divergent part p < gives at most a negli- 
gible contribution to our measured Xt- Also, as 
(3 increases, the physical charge is spread over a 
larger lattice volume and the short range peak of 
n{rmin) becomes more pronounced relative to the 
background exponential, allowing a cleaner filter- 
ing of dislocation-induced charges. Upon remov- 
ing these charges we find that Xt changes by about 
1/3 of the statistical uncertainty for /3 = 6.0 and 
only about 1/10 of the uncertainty for /3 = 6.4. 
We take this as evidence that low-action disloca- 
tions have little effect on our measured suscepti- 
bility, and present results without removing ques- 
tionable charges. The decrease of {Nq)/{Q^) and 
the small effect of questionable charges on mea- 
sured Xt are also observed for other values of the 
RDC cutoff. A detailed analysis of the RDC cut- 
off will be presented in a forthcoming paper [ p^ . 

Since our calculation of xt spans a wide range 
of (3 we can check for scaling by plotting di- 
mensionless ratios of to four other quanti- 
ties versus (3. Our computed values for x^^a 
are 0.114(8) at /3 = 6.0 (16^ lattice), 0.110(14) 
at /5 = 6.0 (243 lattice), 0.090(6) at /? = 6.2, 
and 0.072(5) at /3 = 6.4. The height of each 
set of points is irrelevant since our concern is 
whether the ratios are independent of [3. Scaling 
with nip, /jr, and Kiatt derived from the two-loop 



Table 2 



Lattice Observables (Lattice Units) 



13 a^/'^a nripa f^^a 


6.0 0.220(2)[ 
6.2 0.158(1)1' 
6.4 0.119(2)1' 


1] 0.333(8)[; 
i] 0.277(25)1 : 
i] 0.211(7)1' 


]\ 0.056(15)[31 
^1 0.044(7)1 2] 
i] 0.040(8)1 4] 



formula is apparent, but a possible slope exists 
in Xt'^ l^f^- We present physical values derived 
from nip = 770 MeV and ^ = 440 MeV in Ta- 
ble |l|. Although the issue of scaling still must 
be resolved, we see that these results are consis- 
tent with earlier results by Gockeler et al.\ |l2[ , 
computed on much smaller lattices. Consistency 
between the smaller and larger lattices at /3 = 6.0 
indicates that our lattices are large enough to con- 
trol finite size effects. We conclude that the ge- 
ometric method yields Xt about a factor of four 
higher than both the Witten-Veneziano predic- 
tion and the xt ~ (180 MeV)"* obtained by cool- 
ing [ ^, and that the superficially divergent part 
of Eq. (||) is unlikely to cause this discrepancy. 
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